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Abstract 

We introduce a measure of quantum correlations in the A^-qubit quantum system 

which is invariant with respect to the SU{2^) group of transformations of this system. 

This measure is a modification of the quantum discord introduced earher and is referred 

to as the unitary or 5'C/(2 )-invariant discord. Since the evolution of a quantum system 

is equivalent to the proper unitary transformation, the introduced measure is an integral 

of motion and is completely defined by eigenvalues of the density matrix. As far as the 

calculation of the unitary invariant discord is rather complicated computational problem, 

^ ' we propose its modification which may be found in a simpler way. The case N = 2 

^>mJ . is considered in details. In particular, it is shown that the modified S'C/(4)-invariant 

lO I discord reaches the maximum value for a pure state. A geometric measure of the unitary 

invariant discord of an A^-qubit state is introduced and a simple formula for this measure 

is derived, which allows one to consider this measure as a witness of quantum correlations. 

The relation of the unitary invariant discord with the quantum state transfer along the 

spin chain is considered. We also compare the modified 5t/(4)-invariant discord with the 

^ ■ geometric measure of 5'C/(4)-invariant discord of the two-qubit systems in the thermal 

3 ! equilibrium states governed by the different Hamiltonians. 

(v-i 1 Introduction 

>. 

^ . The development of the quantum information technology stimulates a deep study of the prop- 
t^^ ! erties of quantum correlations inherent in a quantum system. In particular, there is a problem 
"^ I of identification of those quantum correlations which are responsible for the advantages of the 
t^^ ' quantum computations in comparison with the classical ones. The entanglement [H El El HI [5] , 
which was originally taken as a measure of such correlations, seamed out to not cover all of 
them. As a consequence, there are quantum systems without entanglement, which, nevertheless, 
either reveal a quantum nonlocality P, El [8] or speed-up certain calculations in comparison 
K> [ with the classical analogues P, [IDl [HI [H [13]. Such observations cause a new stimulus for 
^ ; study those quantum correlations which are not captured by entanglemet. Thus, the con- 
cept of quantum discord is intensively developing diring last years [TH [151 iMj • Originally the 
quantum discord was introduced to characterize the impact of the classical measurements on a 
quantum system with the purpose to get the maximal information about this system with the 
minimal influence on it [IT]. At flrst glance, the quantum discord seems to cover all quantum 
correlations. However, it was shown that the quantum discord may be either bigger or smaller 
then the entanglement [T5l [T8l [T9] . Thus, we may state that the quantum discord involves dif- 
ferent quantum correlations then entanglement does, in general, which causes doubts whether 
the discord captures all correlations. Moreover, the discord is not symmetrical with respect 
to the subsystem chosen for the projective measurements [H], which stimulates the futher 
study of this ambiguity [20] and suggests to introduce the discord under two-side projective 
measurements [21] (a symmetrical discord). In addition, both the entanglement and discord of 

1 



O 



a 



bipartite A^-qubit system are not invariant with respect to SU(2'^) transformations. On the 
one hand, this allows the evolution of discord which is important if we are interested in the 
strong quantum correlations between the two chosen subsystems. Owing to the discord evo- 
lution, two originally uncorrelated subsystems may evolve to strongly correlated ones. On the 
other hand, SU{2^) transformation means nothing but a transformation to a new basis, which 
may be written in terms of the states of some new "virtual" particles. The discord between the 
new virtual particles may differ from the discord between the original ones. However, it is not 
clear which discord (either between the original particles or between the virtual ones) prevails 
yielding advantages of the quantum devices in comparison with their classical analogues. Of 
course, if one has to effect the first subsystem by means of the second one, then namely the 
quantum correlations between these two subsystems are important so that the usual discord 
may be a proper measure of these correlations. But if one considers a quantum process inside 
of the quantum system, then the quantum correlations between the above mentioned "virtual" 
particles might be more important resulting to advantages of quantum devices. 

Another motivation for the modification of the measure of quantum correlations is a re- 
vealing such quantum systems which have no quantumness. Let us remember that not any 
quantum system exhibits quantum correlations measured in terms of either entanglement or 
discord. Moreover, the family of quantum systems possessing quantum correlations is still not 
completely characterized. Thus, considering the entanglement as a proper measure of quantum 
correlations in a bipartite system (a system composed by two subsystems A and B) [I] , it was 
demonstrated that separable states, i.e. states representable in the form 

i i 

have zero entanglement and, consequently, may not be considered as candidates for implemen- 
tation in quantum devices. Here p^ and p^ represent the density matrices in the Hilbert 
spaces of respectively subsystems A and B. However, the discord as an alternative measure of 
quantum correlations [13] may be nonzero even for separable states ([1]). It was shown [TH [2T] 
that the discord under one-side projective measurements on subsystem B is zero if pa,b is 
representable in the following form 

PA,B = ^PiP'A\bi){hil ^Pi = l, (2) 

i i 

where p^ , 2 = 1,2,..., are density matrices in the Hilbert space of the subsystem A and vectors 
\bi), i = 1,2... represent an orthonormal basis in the Hilbert space of subsystem B. For the 
quantum discord under two-side projective measurements [21], it is shown that this discord is 
zero only if a state is representable in the form 

PA,B = ^Pij\ai)\bj){a^\{bj\, Pij > 0, J^pi^ = 1, (3) 

where vectors \ai), i = 1,2, ... , and \bj), j = 1,2,..., are some orthonormal bases in the Hilbert 
spaces of respectively subsystems A and B. Both eqs.([2]) and ([3]) deffer from eq.([T]) in general. 
Eq.([2]) means that only states diagonalizable by local transformations (i.e. transformations 
in the Hilbert spaces of subsystems A and B) have zero symmetrical quantum discord. Zero 
quantum discord means that all mutual information encoded in the system AU B may be 
revealed by classical measurements using a proper complete set of projective measurements. 



However, the fact that not any complete set of projective measurements may be used to reveal 
all mutual information may be considered as a quantumness of the system. In other words 
one can propose that the system has no quantumness if only any complete set of projective 
measurements exhibits all mutual information encoded in the system AU B. 

All these arguments suggest us to modify (or generalize) the definition of the quantum 
discord. Namely, let us introduce a measure of quantum correlations which takes into account 
the quantum correlations between all possible virtual spin-1/2 particles in the iV-spin system. 
We refer to this measure as the unitary invariant discord and assume that namely this measure of 
quantum correlations is responsible for the advantages of quantum computations in comparison 
with the classical ones. The family of states with zero unitary invariant discord is revealed in 
Sec J4.1[ It is shown that the new measure may be zero only for quantum states whose density 
matrix is proportional to the unit one. 

It is seemed out that the calculation of unitary invariant discord is a complicated computa- 
tional problem involving both the multidimensional optimization and the multiple integration, 
see Sec J2.1[ Therefore we introduce a modification of the unitary invariant discord which re- 
quires only the multiple integration without optimization. In addition, a geometric measure 
[22| [23| [211 [25] of the unitary invariant discord is proposed, which seemed out to coinside with 
the geometric measure of modified unitary invariant discord. The feature of the introduced 
geometric measure is that its calculation does not involve neither optimization nor multiple 
integration. It has a simple analytical representation in terms of eigenvalues of the density 
matrix. For this reason, the geometric measure may be taken as a witness of quantum cor- 
relations in a system. The simplicity of the geometric measure of unitary invariant discord is 
especially important because the calculation of usual descord [H] is a complicated procedure 
which has been carried out only for some particular states, see for instance [151 [El [2T| . 

Hereafter, a system of N spin-1/2 particles is referred to as A^-qubit system. This paper is 
organized as follows. The unitary invariant discord for a two spin-1/2 particle state {SU{A)- 
invariant discord) is introduced and studied in SecJH A modification of the S'f/(4)-invariant 
discord is suggested as a relatively simple method to estimate this measure of quantum correla- 
tions in Sec J2.1.11 It is shown that this modification reaches its maximal value for a pure state. 
The concept of unitary invariant discord is generalized for a system of A^ > 2 spin-1/2 particles 
(5'f/(2^)-invariant discord) in SecJHl A geometric measure of the unitary invariant discord is 
proposed in Sec|H Comparison of the modified S'?7(4)-invariant discord with the normalized 
geometric measure of the S'?7(4)-invariant discord is given in Secj5]for the thermal equilibrium 
states with different Hamiltonians. The basic results are collected in Secj6l 

2 A system of two spin-1/2 particles 

Let us consider a system of two spin-1/2 particles A and B (a two-qubit system) in the standard 
multiplicative basis associated with the above two particles: 

|«i) = |00), |a2) = |01), |«3) = |10), |a4) = |ll). (4) 

Here 1 and mean the spin directed along and opposite the fixed 2;-axis respectively. Let us 
consider an arbitrary special unitary transformation U{(f) G 5't/(4), where the set of indepen- 
dent parameters ip = {ipi, . . . , ipi^) parametrizes the group S'[/(4) with f/(0) = J4 (hereafter /„ 



is the n X n identity matrix). This transformation allows us to introduce another basis 
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|A((^)) = J]t/*.(^)|a,), |A(0)) = |a.), 2 = 1,. ..,4. (5) 

i=i 

The basis |/3j (</?)) may be treated as the multiplicative basis of two virtual spin- 1/2 particles 
A'^ and i?*^. In other words, we may write vectors \l3i{ip)) in the following form similar to the 
representation (jl]): 

|/3i((^)) = |00)^, I/32M) = |01)^, I/33M) = |10)^, |/34M) = |11)^. (6) 

Now the density matrix p for the above system may be written using either the basis \ai) or 
the basis \(3i{ip)), z = 1, . . . , 4, as follows: 

4 4 

p = E p^M\^^){^^ = E p..MiAM)(/3,Mi, iA(o)) = i«.). (7) 

Introduce the matrices p{ip) = {pij{f)} and U{ip) = {Uij{ip)}. The transformation betwen 
bases |q!j) and |/9j(v?)) given by eq.([5]) yields the evident relation: 

p{cp) = f/Mp(0)f/+(^). (8) 

Here p(0) and p{(p) are the density matrices written in the bases |aj) and \f3i{(p)) respectively 
and, consequently, these matrices describe the systems of two particles AU B and A''^ U B'^ 
respectively. Of course, the discord between the particles A and B differs from the discord 
between the (virtual) particles A*^ and -B*^. However, as it is mentioned in the Introduction, both 
discords measure the quantum correlations in the same two-qubit system and it is impossible to 
predict which discord is really responsible for the advantages of quantum devices in comparison 
with their classical counterparts. This is a motivation to introduce a measure of quantum 
correlations which takes into account correlations between all virtual particles A*^ and B"^. 

Remark that relation ([H]) may be viewed as a transformation of the matrix p(0) itself rather 
then the transformation of the basis. This means that the discord between two virtual particles 
A'^ and B'^ is equivalent to the discord between the original particles A and B whose state is 
transformed by the unitary transformation U{ip): U {(p) p{0)U~^ {(p) . This evident remark will 
be used for the calculation of the discord between the virtual particles A'^ and B'^ hereafter. 

2.1 The SU(4) invariant discord and its modification 

The quantum discord as a measure of quantum correlations in a bipartite system has been 
introduced in ref.[Tl]. This definition is based on the projective measurements over one of the 
subsystems. Since there is an asymmetry with respect to the subsystem chosen for the projective 
measurements p3l [20] , a two-side measurement generalization of the bipartite discord has been 
proposed in ref. [21] where the above asymmetry disappears. We use the discord under the 
two-side projective measurements hereafter. 

Let us recall the formula defining the discord for the system consisting of two spin- 1/2 
particles A'^ and B"^ with the density matrix p{ip) [HI [21] : 

QABipiv)) = Hpiv)) - sup /(p((p))- (9) 

{n} 



Here 

I{p{v)) = Sip^i^)) + Sip^i^)) - Sipi^)), (10) 

and p{(p) (corresponding to the two-side projective measurements [21]) is defined as follows: 
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p{^) = J2 %(x)p(v.)n,,(x), (11) 

where Ilij{>i) have the following form: 

nij(x) = Ui{>CA)®Ilj{>CB), >c={xa,Xb), i,j = 1,2 (12) 

and 

n,(x^) = V{ka)MO)V+{xa), M^b) = V{xB)MO)V^i^B), (13) 

v{>€a), v{hb) e su{2), ni(o) = |o)(o|, n2(o) = |i)(i|. 

The vector parameters ka = {xAi, i = 1,2,3} and xb = {>iBi, i = 1,2,3} represent two 
different parametrizations of the group SU{2). Of course, these parameters depend on </?, i.e. 
>^Ai = >^Ai{f) and xsi = XBi{'^)- The optimization in eq.([9]) yields the proper values for the 
parameters xa,Bi which, in turn, depend on ip. quantumness This optimization provides the 
invariance of the discord with respect to the local SU{2) transformations of particles A and B. 
Next, in order to take into account the discords between all possible virtual particles A"^ and 
B"^ we average qAB{p{f)) over ip. However, before proceed to calculations, let us recall that 
there is a freedom in the choice of the initial basis (i.e. the basis corresponding to (f = 0) and 
the result depends on this initial basis. To remove this ambiguity we suggest to take the basis 
of eigenvectors of p arranged in the order of decreasing eigenvalues as the original basis for the 
subsequent calculations. The density matrix p is diagonal in this basis: 

3 

p(0)=A2 = diag{Ai,A2,A3,A4}, A4 = l-^Ai, (14) 

i=l 

< A, < 1, Ai > A2 > A3 > A4. 
Thus, Ai is the maximal eigenvalue hereafter. Using eq. (ll4l) . we rewrite eq.([8]) as 

p(^) = Uicp)A,U+{cp). (15) 

Now the averaged value of qABipi'^)) can be calculated as follows: 

gAB(A2) = f dniip)qAB{p{^)) = (16) 

dn{^)I{p{ip))-snp /"dfi(<^)/(^H,,(x(<^))p((^)H,,(x(<^))), 



u 



where ^2 is a measure in the space of parameters (fi, J dQ{(f) = 1 and Hjj are defined in 
eq. (TT2|) . Remark, that the averaged discord qab may be simply calculated for a pure state 
when qAB^pi,'^)) = S{p^{(p)) = S{p^{(p)) [12]. In this case, there is no optimization in eq. (fT6ll 
and one obtains qAB{^2^^^) ~ 0.42, where 

Ar'^ = diag{l, 0,0,0}. (17) 

Now we represent a list of important properties of qAB{.-^2)- 



1. Due to the averaging, the quantity qAB{^2) is completely defined by the eigenvalues of 
the matrix p. 

2. As a consequence of the previous property, qAB{^2) is invariant with respect to SU{A) 
transformations of the whole system. 

3. Since discord qAB > 0, one has that ^^5(^2) > by its definition f lT6|) . 

4. The averaged discord qab is an integral of motion, i.e. it does not evolve. In fact, any 
evolution is described by the Liouville equation iptit) = [H{t),p{t)] {h = 1) with some 
initial condition p\t=o = Po and the Hamiltonian H. Thus, the evolution of the density 
matrix p may be written as p{t) = U{t)pQU~^{t), U{t) G SU{4). This means that the 
density matrix p{t) has the same eigenvalues as po and consequently the unitary invariant 
discord for p(t) equals to that for po- 

A consequence of the property 2 is that qab is invariant with respect to the local SU{2) 
transformations as well, similar to the usual discord. 

The averaged discord qAB{-^2) can be normalized as follows: 

<5ab(A2) = , .pures ) (18) 

qAB[J^2 ) 

SO that Q abI-^^!^^'^) = 1- The quantity Qab{^2) is referred to as the unitary invariant discord 
hereafter. It possesses the same properties 1 - 4 as the averaged discord ^^5(^2)- 

2.1.1 Modified unitary invariant discord 

The calculation of qAB{-^2) is a complicated computational problem because the second integral 
in eq.( fT6l) involves an optimization over >i{(p). Thus, it is reasonable to introduce another 
quantity which would be simpler for evaluation and gives a reasonable estimation of 5^5(^2). 
The upper bound of (Jab seems to be such quantity. This upper bound may be readily written 
due to the obvious inequality 

sup /"rfi](v9)/(5^n,,(x(^)p(^)n,,(x(^))) > /"df](^)/(5^n,,(o)p(^)n,,(o)),(i9) 

ij ij 

which suggests us to define the upper bound of (Jab as follows: 

Qab{^2) = [ dn{^)i{p{^)) - I dn{^)i{p\^)), p\^) = ^n,,(o)p(^)n,,(o), (20) 

where p'^(v5) is the diagonal part of p{ip). Similar to g^B, its upper bound q^^ depends only on 
the eigenvalues of the matrix p(0). 

One has to emphasize that the upper bound g^'^(A2) may be viewed as a difference between 
the total (i.e. quantum and classical) mutual information encoded in the two-qubit system 
and its classical contribution encoded in the diagonal part of the matrix p(</3) for all (p. From 
this point of view, the value g^^(A2) may be considered as an alternative measure of quantum 
correlations. Then it is reasonable to normalize q^^ as follows 

QaB\^2) = -up f i.pure\ i (21) 

9ab V^2 ) 



so that QABi,^^"^^) = Qab{-^2^^'^) = 1- By construction, the value Qab is not an upper 
bound of the unitary invariant discord Qab in general. We refer to Qab{-^2) as the modified 
unitary invariant discord. Emphasize that both the normalization ^^^^(Ag"'^) in eq.( !T8|) and 
the normalization ^^^(Ag*"^^) in eg. (1211) are the universal constants, i.e. they remain the same 
for all two-qubit systems. 

It will be shown below in this section that ^^^(Ag*"^^) ~ 0.69 and Qab{^2) < 1- In 
addition, it is obvious, that the modified unitary invariant discord possesses the same list of 
properties 1 — 4 as the unitary invariant discord. 

Finally remark that there is only one possible bipartite decomposition of the system of two 
spin-1/2 particles, so that we may write 

g(A2) = Qab{A2), Q(A2) = Qab{A2). (22) 

Modified S't/(4)-invariant discord as a function of density matrix eigenvalues. We 

investigate the A2-dependence of the modified unitary invariant discord Q(A2) in this paragraph. 
This will give us a clue to reveal systems with strong quantum correlations. For this purpose, 
let us use the parametrization of an arbitrary U{(f) G SU{4) proposed in ref. |26j : 

JJtm) = g*73</'lg*72¥'2gi73</'3gn5V4gi73V5gni0V6gn3</'7 ^ /'23~) 

gn2</'8gi73¥'9gi75</'10gi73'/'llgn2</'12gn3</'13gi78V14g«715'/'15 

where 4x4 matrices ji are given in the Appendix, eqs. fl52l) . We parametrize A2 by three 
parameters Aj, i = 1,2, 3, as it is shown in eq. flT^ . Eq. fll4l) represents the most general structure 
of matrix A2 with the eigenvalues Aj on the diagonal arranged in the decreasing order. Note 
that not all parameters (pi appear in the matrix p{(p) = U{(p)A2U~^{(p). In fact, 73, 73 and 715 
are diagonal matrices (see eqs. flS^ ) so that they commute with each other and with A2. Thus, 
the formula for p{ip) reads: 

piv) = U{ip)K2U^{ip) = U{ip)K2U^{ip), (24) 

IJ(in) = g«73'/'lg«72¥'2g«73'/'3g«75¥'4gi73'/'5g«7l0¥'6g«73'/'7g«72l/58g«73'/'9g«75¥'10g«73'/'llg«72'/'12 

which possesses 12 parameters cpi, i = 1, . . . , 12. Here ranges of parameters ipt are following 

m 

< fl, f3, ^5, fr, f9, V'll < TT, < V92, (fA, V^G, ^8, V^IO, V^12 < ", (25) 

2^ 
so that dQ(ip) = -—rdipi . . . d(pi2. 

To analyze the A2-dependence of the modified unitary invariant discord Q, we, first of all, 
consider the case of at most two nonzero eigenvalues A^: A2 = 1 — Ai, A3 = A4 = 0. Thus, 
there is one independent parameter Ai in this case. The graph of the function Q(A2) = Q(Ai) 
is shown in FigJU^a). Here Q > 0.72, i.e. the modified S'?7(4)-invariant discord may not be less 
then ^ 0.72 if there are two zero eigenvalues. 

Next, consider the case of at most three nonzero A^: A3 = 1 — Ai — A2, A4 = 0. Thus, 
there are two independent parameters Ai and A2 in this case. The graphical representation 
of the function Q(A2) = Q(Ai,A2) is shown in FigHt^fc). Here Q > 0.30, i.e. the modified 
S'f/(4)-invariant discord may not be less then ^ 0.30 if there is one zero eigenvalue. 
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Figure 1: The modified S'f/(4)-invariant discord as a function of tlie maximal eigenvalue Ai 
for different numbers of nonzero eigenvalues; (a) at most two nonzero eigenvalues Ai and A2: 
A2 = 1 — Ai, A3 = A4 = 0; Q > 0.72; (fe) at most three nonzero eigenvalues Aj, i = 1,2,3: 
A3 = 1 — Ai — A2, A4 = 0; Q > 0.30; (c) all eigenvalues may be nonzero: A4 = 1 — Ai — A2 — A3; 
Q>0 



Finally, the most general case of all nonzero eigenvalues is represented in Fig{T]^c). In this 
case Q(A2) = Q(Ai, A2, A3) > 0. Note that the modified S't/(4)-invariant discord (as well as the 
5't/(4)-invariant discord itself) equals zero only in the case A2 = jh, see Sec J4.1l and FigHJ^c), 
the curve with A2 = A3 = 0.25. As opposite, there is a family of states where the usual discord 
is zero [H]. FigH] demonstrates that the modified 5't/(4)- invariant discord does not exceed unit 
(Q < 1) with ^(Aj"^^) = 1. In general, we conclude that in order to maximize the modified 
S'f/(4)-invariant discord, one has to build the density matrix with as many as possible vanishing 
eigenvalues. In particular, the condition (1 — Ai) ^ 1 may serve as a sufficient condition for 
the large modified S'f/(4)-invariant discord. 

3 Bipartite unitary invariant discord in an A^-qubit sys- 
tem 

Let us introduce the unitary invariant discord for the A^-qubit system, N > 2 (or for the 
system of A^ > 2 spin- 1/2 particles). Similar to the two-qubit system, we take the basis of the 
eigenvectors of the density matrix arranged in the decreasing order of its eigenvalues as the 
original basis (i.e. the basis corresponding to c/? = 0). The density matrix is diagonal in this 
basis: 

2^-1 
p(0) = A^ = diag{Ai,...,A2^}, A2iv = l-^A„ (26) 



i=l 



0<A, <1, Ai > A2> ■■■> A2iv, 



which is the A^-qubit analogy of eq.( 1T^ . Next, using this basis, we consider any two subsystems 
A and B with A^^ and Nb spin-1/2 particles respectively, Na + Nb = N. Now we may introduce 
the 5'f/(2^)-invariant bipartite discord Qab and its modification Qab just using formulas (^ 
^T]) with the replacement A2 — )■ A^r and with the vector parameter ip having the appropriate 
dimension: (p = {ipi, . . . , (Pmn)^ ^n = 2^^ — 1. The two-side projectors Ilij will be generalized 



IN n 

for the multiqubit case as follows: 



n,,(x) = nf(xA)®nf(xB), x=(xa,xb), (27) 

where >iA,B are vector parameters ka = {x^i, . . . , x^j-g^^A-i)}, ><b = {x^si, ■ ■ ■ ,Xb{2'^'^b-i)} 
and 

nf(x^) = \/^(x^)nf(0)(O + (xA), F^(x^) G 5f/(2^-^), z = l,...,2^-, (28) 
Iif{xB) = V\xB)IifmV^Y{>CB), y^(xB) G ^f/(2^-), 2 = 1,..., 2^-. 

Each of the projection operators 11^ ' (0), z = 1, . . . , 2^-^'^, picks up the zth diagonal element 
of the proper subsystem: 

nf(0) = |n,,i...ni,^^)(ni,i...n,,;v^|, 2 = 1,...,2^^, (29) 

nf (0) = |nj,i . . . nj,7v^)(raj,i . . . nj,7v^|, j = 1, . . . , 2^^, 

where n, .,■ equal either or 1 and |nj^i . . . Ui^N^) and |nj i . . . Uj^Ng) are bases in the Hilbert sub- 
spaces corresponding to the particles A and B. The multiqubit generalization of the eigenvalue 
matrix A^^^'^ (see eq.( !T7|) ) is following: 

Ar'^ = diag(l,0,...,0,0). (30) 



In eqs.( p6|) and (!30|) . the eigenvalue Ai corresponds to the zero eigenvectors | ... 0) and | ^. . 0). 

It is obvious that the unitary invariant discord Qab (as well as the modified unitary invari- 
ant discord Qab) written for any particular choice of A and B does not capture all possible 
bipartite quantum correlations in the system of N particles. To take into account all possible 
bipartite quantum correlations we must consider the invariant discords of all possible bipartite 
decompositions. Namely, we suggest the following definitions of the bipartite unitary invariant 
discord and the modified bipartite unitary invariant discord of the A^-qubit quantum system: 

Q{An) = min QAsi^N), Q{^n) = min Qab{^n), (31) 

A,B A,B 

where minimization is taken over all possible bipartite decompositions of the quantum system 
AUB. Considering the system of two spin-1/2 particles we see that the minimization disappears 
from eqs. fl3T]) so that we obtain eqs. fl22l) . Finally, let us note that the definition of the (modified) 
bipartite unitary invariant discord given by eqs. fl3Tl) is similar to the characteristics of the 
multipartite entanglement introduced in ref. [27]: both use the minimization over all possible 
bipartite decompositions of a quantum system. 

Discords Q{A]y) and Q^A^) have the same four properties as respectively Q(A2) and Q(A2) 
(see Sec J2.ip with replacements 

As -^ An, SU(4) -^ SU{2^). (32) 



3.1 The bipartite unitary invariant discord of the reduced A^-qubit 
system 

Now we turn to the discord between subsystems A and B consisting of A^^ and Nb particles 
with Na + Nb < N. As usual, the first step is the proper reduction of the original density 
matrix with respect to the subsystem C consisting of Nc = N — Na — Nb particles. However, 
considering the unitary invariant discord, we are not forced to perform the reduction in the 
system of original physical particles unlike the usual discord. Instead of physical particles, we 
are free to choose a proper set of the virtual particles representing the given quantum system 
and perform the reduction inside of this set. As it was written in the Introduction, the choice of 
the virtual particles is equivalent to the choice of the proper transformation U{ipi) G SU{2'^), 
which transforms the original basis (which is the basis of the eigenvectors of p arranged in the 
decreasing order of eigenvalues) into the multiplicative basis of the appropriate virtual particles, 
where the density matrix reads: p{(p^) = U {(p^) AnU'^ {(p^) . Now we reduce p{(p^) with respect to 
the appropriate subsystem C and calculate the unitary invariant discord Qc using formulas of 
Sec J3] with the reduced density matrix having dimension 2^^+^^. Here Qc means the bipartite 
unitary invariant discord of the quantum system reduced with respect to the subsystem C. 

There is a natural question whether the calculation of the discord in the reduced quantum 
system is a valuable procedure. In fact, the unitary invariant discord Q given by the first eq. fl3ip 
takes into account all bipartite correlations in a quantum system. Nevertheless, the unitary 
invariant discord in a reduced quantum system may be useful, for instance, in the following 
two cases. 

1. We are interested in the discord in the subsystem of A^^ + A'^^ < N original particles. In 
this case we consider the density matrix of quantum system in the multiplicative basis 
corresponding to the original particles and reduce it with respect to the subsystem C. 
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2. As far as the reduced density matrix has lower dimension, the calculation of the bipartite 
unitary invariant discord for this matrix is simpler. This fact may be used to obtain an 
information about the overall quantum correlations in a large quantum system having the 
bipartite unitary invariant discords of the appropriate reduced subsystems, as it is done 
in the Sec J3.1.11 

3.1.1 bi-particle unitary invariant discord of an A^-qubit spin-1/2 system in a pure 
state 

We consider the case when the calculation of the unitary invariant quantum discord in the 
reduced quantum system is helphull to demonstrate the presence of strong bipartite quantum 
correlations. Let us consider the single quantum state transfer along a spin chain [281 EHl EDI 
[3T1 l32] and turn to the problem of relation between the single quantum state transfer and 
the quantum correlations. There are articles where the relation between the state transfer 
probability and the entanglement is studied [331 EH ES]. It is demonstrated that the perfect 
end-to-end state transfer probability and strong entanglement are not necessary related with 
each other. It has also been shown recently that one can arrange either the high probability 
quantum state transfer or the strong entanglement between any two particles in the spin-1/2 
chain of A^ nodes governed by the Hxyz Hamiltonian with the inhomogeneous magnetic field 
[35] (similar numerical experiments with discord have not been evaluated yet). This fact may be 
taken as an evidence of the quantum correlations which supplement the quantum state transfers 
in the above systems. Of course, these correlations may be caught by the unitary invariant 
discord (IHTl) . However, calculations are very complicated because of the numerous parameters 
in SU{2^). For this reasong we use the advantage of the reduced density matrix. Namely, 
we show that an arbitrary system of A^ spin-1/2 particles in a pure state may be viewed as a 
system of A^ virtual spin-1/2 particles with 

'^ij;rest ^ij;rest ^ \^^ ) 

for any two virtual particles i and j. Here Qij-rest and Qij-rest niean the S'?7(4)-invariant discord 
and the modified S'f/(4)-invariant discord between virtual particles i and j in a system of A" 
virtual spin-1/2 particles reduced with respect to all other virtual particles. The above system 
of virtual particles corresponds to the basis of ordered eigenvectors of the density matrix, i.e. 
p(0) = A^"'^^ see eq.dSOj). 

First of all remember that the only nonzero eigenvalue Ai = 1 corresponds to the zero 
eigenvectors | ... 0) and | ... 0), as it is pointed after eq. (l30l) . It is simple to show that the 

reduced density matrix with respect to any (A^ — 2) virtual particles reads 

p-'^ = A2 = diag(l, 0,0,0), (34) 

which is the density matrix of a pure state (fTT]) . Thus Qij-rest{-^2) = Qij;rest{-^2) = 1 for any 
two particles i and j, and eqs.(l33l) are valid. This allows us to conclude that the bipartite 
quantum correlations in the whole spin system are big. 

This conclusion is applicable to the process of a single quantum state transfer along the 
spin-1/2 chain of A^ nodes since the initial state is a pure one in this process. Thus we conclude 
that a single quantum state transfer along a spin chain is directly related with strong quantum 
correlations in this spin chain. 
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4 Geometric measure of a bipartite unitary invariant 
discord 

A geometric measure of either quantum entanglement or quantum discord has been introduced 
as a characteristics showing how far is a given state from that having zero entanglement or 
discord. It is defined by the formula [22| [23| l2^ [25] 

lABi^N) = J dn{^) inf I Ip(^) -x\\\ (35) 

where x is the complete set of density matrices having zero measure of quantum correlations 
(either entanglement or discord). Here ||p|P = Tr(p^). In our case, x is the complete set of 
matrices having such ordered matrices of eigenvalues A^ that Q{A^) = 0. Thus, first of all, 
one has to define the appropriate set of matrices x- 

4.1 Quantum state with zero bipartite unitary invariant discord 

Let us show that the bipartite unitary invariant discord introduced by the first formula (!3T!) is 
zero only if all eigenvalues of the density matrix equal each other, i.e. 

X = ^N = ^V- (36) 

First of all, remember that the usual discord is zero for the density matrix p{(po) which may be 
diagonalized by the local unitary transformations V"^ G SU{2^^) and V^ G SU{2^'') p^[2T] : 

p(^o) = V^ixAiifo)) ® V''{Mvo))^n{V^V{M^o)) ® (\/^)+(xb(^o)), (37) 

where A^v is a diagonal matrix of eigenvalues of p{(po)- Note that the eigenvalues in A^ are 
not necessary ordered, but they may be ordered by the appropriate orthogonal transformation 
O e S0{2^y. 

An = O+AnO. (38) 

Here A^v is an ordered matrix of eigenvalues ( 126|) . Regarding the unitary invariant discord, it 
may be zero only if eq.( l37|) is valid for any (fQ, which is a consequence of the nonnegativity of 
discord. Let us find the set of such matrices A^ that satisfy this condition. As usual, we start 
with the basis of eigenvectors of p which are arranged in the decreasing order of eigenvalues, i.e. 
p(0) = Atv and p{(p) = U{(p)ANU~^{(f). Then the unitary invariant discord Q{An) is zero if, for 
any ip, there is such local transformation V^{xa{'^)) G SU{2^-'^) and V^^xsi^p)) G S'f/(2^^) 
that 

V^{>ca{^)) ® y^(xB(<^))f/(^)A;vf/+(<^)(V^^) + (xA(<^)) ® {V'yixBi^)) = An, (39) 

where A^v consists of the eigenvalues of p(0) and, consequently, there is an orthogonal trans- 
formation O G 5*0(2^) such that A^v = O^AnO. As a consequence, eq.(l39l) may be written in 
the following form: 

[V^ixAiif)) ® V''{HB{ip))U{ip)0+, An] = 0. (40) 
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Since [/(</?) E SU{2^), then the matrix V^^(xa(v5)) ® V^{HBi^))U{(^)0+ may not be diagonal 
for any if. Thus, in order to satisfy condition (140|) . one has to assume that A^v is proportional 
to the identity matrix, i.e. 

AtV = ^/2iv (41) 

independently on the particular choice of subsystems A and B. This means that the state 
with zero unitary invariant discord is the state having all equal eigenvalues, i.e. eq. (l36|l is 
valid. Then the manifold of matrices x (see eq.( l35|l ) with fixed A^ consists of the single element 
X = ^l2^- This conclusion allows one to simplify eq. (l5H]) . 

In addition, one can show that the modified unitary invariant discord (see the second formula 
fl3T]) ) is zero only for the density matrix fl36|) similar to the unitary invariant discord. In fact, 
since q^^ > qab, then g^^ may be zero only if qab = 0. Thus, one has to show that q^B is 
zero for the density matrix ( 136|) . The later is true because, if p(0) = x? then p{ip) = x as well, 
which means that p'^{ip) = pif)- Then eq. (l20l) yields q^Bix) = 0- Consequently, Q (defined by 
the second equation ( l22i) ) is zero only for the matrix x given by eq.( l36l) . A direct consequence 
of this conclusion is that the geometric measure of modified unitary invariant discord coinsides 
with the geometric measure of unitary invariant discord. 

4.2 Explicite formula for the geometric measure 

In accordance with Sec J4.lt eq. flS^ gets the following form: 

9ab(A^) = I dQiip)\\p{v) - ^Vir = (42) 

I dQiip)\\Uiip)iA^ - ^V)f/+(c^)||2 = ||A;v - ^VIT- 
Thus 

qU^n) = 5^(A. - ^)^ = E A? - ^, (43) 

where we take into account that ^^ Aj = 1. We see that the geometric measure of the unitary 
invariant discord does not require neither optimization no multiple integration. It is obvious 
that this measure takes the maximal value for the pure state, when there is a single nonzero 
eigenvalue Ai = 1: 

g^(Ar^) = 1 - ^. (44) 

In particular, eq.(g§ with A^ = 2 yields g'^(Af"''^) = 3/4 > g(A^"''^) ^ 0.42. It is reasonable to 
normalize the geometric measure as follows 

00(, ^ _ g^(A^) _ 2"E£iA^-1 .4,. 

thus Q^{A]\f) < 1. The quantity Q'^{Aj^f) is referred to as the normalized geometric measure 
of the unitary invariant discord (or normalized geometric measure for the sake of brevity). It 
possesses the same properties as Q{Aj^), see list of four properties in Sec J2.1l with replacements 
(132 p . It is remarkable that formula f l45p remains correct for an arbitrary spin-1/2 quantum 
system. 
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5 Example: modified S'?7(4)-invariant discord and nor- 
malized geometric measure of thermal equilibrium states 

As an example, let us compare the modified S't/(4)-invariant discord with the normalized 
geometric measure of the system of two spin- 1/2 particles in the thermal equilibrium state 
[361 EH [20] 



Z = Ti{e-kT} 



(46) 



with different Hamiltonians H: 

Hueis 

HxYZ 

HxY 



—D[Iixl2,x + h,yl2,y ~ '^h,zl2,z), 
~D[Ii;^l2^x + h,yl2,y)- 



(47) 
(48) 
(49) 



Here D (either positive or negative) is the constant of either dipole-dipole or exchange inter- 
action, T is the temperature, k is the Boltsman constant. In order to calculate the modified 
unitary invariant discord one needs the eigenvalues of Hamiltonians: 

_ D D D 3D 
D D 

A^,, = {0,0,--,-}. 

The normalized geometric measures of the thermal equilibrium states with Hamiltonians Wf\ 
BHj) may be readily calculated using eq. f HSjl with eigenvalues Aj corresponding to the eigenvalues 
of Hamiltonians ( 150|) : 



Q 



G 
Heis 



QxYZ 



oP 



3e/5 + l 



L>>0 



D<0 



Ve^ + 3. 

3e3/3 _ 2e2/3 _ ^^W/2 + 3^/3 _ 4^/3/2 ^ 4 

3(e3/3/2 + 6/5/2 + 2)2 ' 

4e3/3 _ 4e5/3/2 + 3^2/3 _ 4^3/3/2 _ 26^^ + 3 



Q 



G 

XY 



3(2e3/3/2 + e/3 - 
-(sech^(/3/4)-4sech^(/3/4) 



If 
■3), 



VD, 



(51) 



D>Q 



D <0 



where /3 = D/{kT). However, the modified unitary invariant discord itself may not be calculated 
analytically. Results of numerical calculations are represented in Figj2] (thick lines), where the 
normalized geometrical measures (15T1) are shown for comparison (thin lines). 

The different values of the modified unitary invariant discords of the ground states (T = 0) 
are well explained by the Hamiltonian eigenvalues (l50l) . Since the ground state corresponds 
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In 




-i — I — I — I — I — I — I — I — I — I — r 
0,5 1 1,5 

DiTfiensionless temperature, k T/D 



Figure 2: The thermal equhbrium state of two spin-1/2 particles: the comparison of the modi- 
fied S'[/(4)-invariant discord Q with the normalized geometric measure Q^ for different Hamil- 
tonians. Graph corresponding to the Hamiltonian Hxy does not depend on the sign of D (there 
is a single thick and thin dashed lines) 

to the minimal eigenvalue (s), it is a pure state in most of the considered cases such as Hxy 
(VD ^ 0), Hdz {D > 0) and Hueis {D < 0), which correspond to Q = Q^ = 1. The ground 
state has the degeneration degree two in the case of the Hamiltonian H^z with D < 0, when the 
nonzero eigenvalues of the density matrix are Ai = A2 = 1/2 and Q ^ 0.72, Q*^ = 1/3, so that 
Q — Q^ ^ 0.39. The ground state has the degeneration degree three in the case of Hamiltonian 
Hheis with D > 0, when the nonzero eigenvalues of the density matrix are Ai = A2 = A3 = 1/3 
and Q ^ 0.30, Q'^ = 1/9 so that Q - Q^ ^ 0.19. The calculated values of Q of ground states 
are confirmed by FigHJ We see that the difference between Q and Q'-^ is significant for the 
degenerate ground states. 

6 Conclusions 

We propose a S't/(2^)-invariant measure of bipartite quantum correlations in a system of A^ 
spin-1/2 particles. The feature of this measure is that it takes into account not only the quantum 
correlations between the original physical particles but also the quantum correlations between 
all virtual particles, which may also be responsible for the advantages of quantum devices. 
The calculation of the unitary invariant discord seemed out to be a complicated computational 
problem involving both the multiparameter optimization and integration over the multidimen- 
sional parameter space of SU{2^) group. This forces us to introduce the modified unitary 
invariant discord which requires only multidimensional integration without optimization and 
consequently may be calculated simpler. The case A^ = 2 is studied in more details. In particu- 
lar, it is obtained that the maximal value of the modified S't/(4)-invariant discord corresponds 
to a pure state. Moreover, considering a pure state of a multiqubit system we obtain that there 
is such complete set of A^ virtual particles that the unitary invariant discord between any two 
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particles from this set achieves unit. This example is closely related with the single quantum 
state transfer along the spin-1/2 chain as an evidence of strong quantum correlations which 
assist this process. This fact confirms the importance of the quantum correlations between the 
virtual particles in the quantum information processes. 

Next, we introduce the normalized geometric measure of unitary invariant discord which 
coinsides with the normalized geometric measure of modified unitary invariant discord. Unlike 
the unitary invariant discord itself, its geometric measure has simple analitical representation 
in terms of eigenvalues of the density matrix and, as a consequence, neither optimization nor 
multidimensional integration is required for its calculation. 

The unitray invariant discord Q, the modified discord Q and the geometric measure Q'^ 
possess the same list of four properties given in Sec J2. II with replacements ( 132|) . In addition, it 
is shown that Q(A2) < 1 and Q'^{An) < 1. 

The calculation of the unitary invariant discord is a complicated procedure involving both 
optimization and multiple integration. The calculation of the modified unitary invariant discord 
is simpler, because it involves only multiple integration. But both of the above discords require 
numerical calculation. However, the calculation of the geometric measure is quite simple and 
may be performed analytically for the states of system with any number of spin-1/2 particles 
using eq.f l45p . This property of the geometric measure allows one to use it as a witness of 
quantum correlations which is especially important because the calculation of the usual discord 
using the algorithm described in [13] is a complicated procedure and can be done only for 
particular states [151 [HI HI]- Calculating the geometric measure using eq. (l45|) we readily 
answer the question whether the given system possesses quantum correlations. 

The modified unitary invariant discords and the normalized geometric measures of both 
the ground states and the thermal equilibrium states for the systems governed by different 
Hamiltonians are compared. FigJ2] shows that the modified 5't/(4)-invariant discords as well 
as the normalized geometric measures of S't/(4)-invariant discord of the ground states (T = 0) 
decrease with the increase in the degeneration degree of the ground states. 

Altough the unitary invariant discord is introduced for the iV-qubit system, it may be 
straightforwardly generalized for an arbitrary quantum system. The same is valid for the 
geometric measure. 

The author thanks Professor E.B.Fel'dman for usefuU discusions and advices. This work 
is supported by the Program of the Presidium of RAS No. 18 "Development of methods of 
obtaining chemical compounds and creation of new materials" . 
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7 Appendix 



The basis for the Lie algebra of SU{A) [38] 
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